Mon. Not. R. Astron. Soc. 000,[T]l6]{201 1) Printed 20 July 201 1 (MN MgC style file v2.2) 



Swift monitoring of the central X-ray source in RCW 103 
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ABSTRACT 

The X-ray source IE 161348-5055 lies at the centre of the 2-kyr-old supernova remnant 
RCW 103. Owing to its 24-ks modulation, orders-of-magnitude flux variability over a few 
months/years, and lack of an obvious optical counterpart, IE 161348-5055 defies assignment 
to any known class of X-ray sources. Starting from April 2006, Swift observed IE 161348- 
5055 with its X-ray telescope for ^2 ks approximately once per month. During the five years 
covered, the source has remained in a quiescent state, with an average observed flux of ^1.7 x 
10"^^ erg cm^^ s^^ (1-10 keV), ^20 times lower than the historical maximum attained in 
its 1999-2000 outburst. The long time-span of the Swift data allows us to obtain an accurate 
measure of the period of IE 161348-5055 [P = 24 030.42(2) s] and to derive the first upper 



limit on its period derivative (|P|< 1.6x10 ^ss ^ at 3 cr). 

Key words: pulsars: general - stars: neutron - X-rays: individual: IE 161348-5055. 



1 INTRODUCTION 

IE 161348-50 55 (hereafter IE 1613) wa s discovered with the Ein- 
stein satellite jTuohv & Garmirelll98(J) close to the geometrical 
cent re of the young supernova remn ant (SNR) RCW 103 (age ~2 
kvr; ICarter. Dickel & Bomanslll997n . It was proposed as the first 
example of a radio-quiet (possibly owing to an unfavourable radio 
beaming), isolated, cooling neutron star i Tuohy & Garmirelll980l ; 



iTuohy et alJl983l : lGotthelf. Petre & Hwanalll997h" 

At p resent, there is little doubt that IE 1613 indeed is a neu- 
tron star jDe Luca et al.ll2006D and the source is traditionally in- 
cluded in the class of the 'central compact objects' (CCOs; see 
lOe Lucall2008l for a review). CCOs are a small group of young 
and seemingly isolated X-ray-emitting neutron stars (with thermal- 
like spectra), observed close to the centre of non-plerionic SNRs 
and without obvious counterparts in other wavebands. However, 
its peculiar temporal behaviour distinguishes IE 1613 from the 
other CCOs (actually, it singles this source out as a unique ob- 
ject in general). The first peculiarity of the source is its orders-of- 
magnitude X-ray flu x variability on a few months/years time-s cale 
JGotthelf. Petre & Vasisht 1999: Ga rmire et al. 2000; Sanwal e t"ai] 
I2OO2I ; iBecker & AschenbachI I2OO2I) . Moreover, the first Chandra 
observation oflE1613ina low state hinted at a possible periodicity 



at ~6 hours JGarmire et alj2000l) that was not confirmed by subse- 
quent observations of the source in bright states. A long (90 ks) ob- 
servation with XMM-Newton, performed in 2005, caught IE 1613 
in a low state and yielded unambiguous evidence for a strong, 
ne arly sinusoidal modu lation at 6.67 ± 0.03 hours (24.0 ±0.1 
ks; lDe Luca et al.l2006l) . The same periodicity was then recognised 
also in the older data-sets, albeit with a very different pulse shape, 
including two naiTow dips per period. No faster pulsations are seen 
in IE 1613 (iPe Luca et al. 2006). 

Large flux variations, similar to those observed in IE 1613, 
are common among magnetars (e.g. iRea & Espositd uOllh . but 
these pulsars, whose emission is believed to be powered mainly 
by the magnetic field, are characterised by rotational periods 
in the narrow range 2-12 s. On the other hand, CCOs are 
steady sour ces, and their periods -when known- are in the 0.1- 
0.5 s range bavlin et alfcOOd ; lOotthelf. Halpern & SewardlbOOJ ; 
iGotthelf & Halpeml2009l) . If IE 1613 is indeed a magnetar, it must 
have been slowed down by some unusual mechanisms, perhaps by 
a pro peller interaction with a debris disk JDe Luca et alj|2006l ; iLiI 
[20071). A different possibility is that IE 1613 is a peculiar low-mass 
binarylj powered by a double (wind plus disk) accretion onto a 
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^ Deep observations of the field of IE 1613 with the Very Lai'ge Telescope 
and the Hubble Space Telescope showed only two or three faint infrared 
sources (H ~ 22) consistent with the position of IE 1613. If none of them 
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recently-born compact object (in this case the 24-ks signal would 
result from the orbital motion of the system) or hosting a magnetar 
toe Luca et alj|2006l : IPizzolato et aPllOOSl : IShadkamkar & Ghosh! 
I2OO9I) . Both scenarios require nonstandard assumptions about the 
formation and evolution of compact objects in supernova explo- 
sions. 

Here we report on the res ults from a 5-year monitoring of 
IE 1613 with the Swift satellite JGehrels et aLl (2004). This unique 
data set allowed us to obtain a phase-coherent timing solution en- 
compassing also Chandra and XMM-Newton archival observations. 
Thanks to this, we are able to derive an accurate period for IE 1613 
and to set the first upper limits on the period derivative for this puz- 
zling source. 
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Figure 1. Swift/XRT light curve of IE 1613 (1-10 keV). Each observation 
is represented by a point in the plot. A rough conversion of 0. 1 count s ~ ^ ~ 
5 X 10^^'^ erg cm~^ s~^ (1—10 keV, not corrected for absorption) can be 
derived from the spectral analysis (Section[3}. 



2 X-RAY OBSERVATIONS AND DATA REDUCTION 

2.1 Swift data 

The X-Ray Telescope (XRT: iBurrows et alj|2005f ) on-board Swift 
uses a front-illuminated CCD detector sensitive to photons between 
0.2 and 10 keV with an effective area of about 110 cm'^ (at 1.5 
keV) and a field view of 23-arcmin in diameter. Two main read- 
out modes are available: photon counting (PC) and windowed tim- 
ing (WT). PC mode provides two dimensional imaging information 
and a 2.5073-s time resolution; in WT mode only one-dimensional 
imaging is preserved, achieving a time resolution of 1.766 ms. 

Between April 2006 and April 2011, IE 1613 was observed 
by XRT 49 times, for a total net exposure time of 102.8 ks in PC 
modelfl The distribution of the Swift observations can be seen in the 
long-term light curve in Fig.[T] while in Fig.|2]we show the image 
of IE 1613 and RCW 103 resulting from all the XRT data gath- 
ered so far Except for periods in which the source was not visible 
by the XRT because of pointing constraints of the Swift spacecraft, 
approximately one 2-ks observation in imaging mode was collected 
per month. On a few occasions, when IE 1613 showed hints of con- 
sistent flux variations, we requested target-of-opportunity observa- 
tions. For example this happened at the end of October 2010 (see 
Fig. [T] around MJD 55 500), when the source count rate remained 
at a relatively high level of Ri0.09 counts s^^ for a few consecu- 
tive pointings spanning ~5 days. The observations were not time- 
constrained, so the monitoring can be considered a casual sampling 
of the phase of IE 1613. 

The XRT data were uniformly processed with XRTPIPELINE 
(version 12, in the HEASOFT software package version 6.9), fil- 
tered and screened with standard criteria. In order to reduce the con- 
tamination from the SNR, the source counts were energy-selected 
in the 1-10 keV band and extracted within a 10-pixel radius (one 
XRT pixel corresponds to about 2'.' 36). To convert the photon ar- 
rival times to the Solar system barycentre for the timing analy- 
sis, we used t he BA RYCORR task and the Chandra position of 
lOe Luca et al.l 1 120081) . For the spectroscopy, we used the latest 
spectral redistribution matrices in CALDB (20091130), while the 
ancillary response files were generated with XRTMKARF, which ac- 
counts for different extraction regions, vignetting and point-spread 
function corrections. 



is linked to the X-ray source, IE 161 3 is undetected in the near infrared 
down to H > 23 JDe Luca et al.l2008l) . 

^ Also 5.5 ks of WT data were collected during the same pointings. How- 
ever, given the bright SNR in which IE 1613 is embedded (see Fig. [2), we 
did not make use of them in this work. 




Figure 2. Swift/XKT image of IE 1613 and its SNR, RCW 103. All the XRT 
data (as of 201 1 April) were combined, totalling to 102.8 ks exposure time. 
Photon energy is colour-coded; red con'esponds to 0.2-0.9 keV energies, 
green to 0.9-1.7 keV, blue to 1.7-8.0 keV. North is up, east is left. The 
source and background regions considered for the analysis (Section [5) are 
over-plotted. 



2.2 Chandra and XMM-Newton data 

We complemented the 5w(/f/XRT data-set with the few Chandra 
and XMM-Newton observations long-enough to contain a mini- 
mum of two modulation cycles of IE 1613 (Table [TJ. The Chan- 
dra/AClS-S observation (performed on 2002 March 03, when 
the source was rather bright) has been already published in 
[Sanwal et al. (2002), the XMM-Newton/EPlC (2005 August 23-24) 
and ChandrafHRC-S (2007 July 03) ones, both ca rried out while 
IE 16 13 was in a low state, have b een published in lOe Luca et al.l 
( 120060 and in lDe Luca et al.l ( |2008|) . respectively; we refer to these 
papers for more details. The Chondral AClS-l observation of 2010 
June 01 is reported here for the first time. For this work the data 
were processed and analysed with standard procedures, using the 
latest available versions of the Chandra Interactive Analysis of Ob- 
servation software (CIAO, version 4.2) and of the XMM-Newton 
Science Analysis Software (SAS, version 10). 
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Table 1. XMM-Newton and Chandra observations used for this work. 



Instrument 


Obs.ID 


Date" 


Duration'' 






(MJD TBD) 


(ks) 


Chandra/AClS-S 


2759 


52336.489 


50.3 


XMM-Newton/EPlC 


0302390101 


53605.824 


87.5 


ChandrafHRC-S 


7619 


54284.330 


80.2 


Cluuidra/AClS-I 


11823 


55348.604 


62.5 



■^ Mid-point of observation. 

^ Time between first and last event. 




Frequency (Hz) 

Figure 3. Power spectrum computed for the combined Swift/XKT data. The 
peak corresponding to the 24-ks signal is evident. 



3 ANALYSIS AND RESULTS 

We merged the data from the Swift/XRT observations and accu- 
mulated a combined spectrum (a detailed spectral analysis will be 
reported elsewhere). The data were rebinned with a minimum of 
20 counts per energy bin and the background was estimated from 
an annular region centred on IE 1613 (with radii 10 and 20 pixel, 
see Fig. |2j. The spectrum can be fit [xl = 0.82 for 174 de- 
grees of freedom (dof)] by a double-blackbody corrected for the 
interstellar absorption (see lOe Luca et al.l 120061) . The best-fit pa- 
rameters are blackbody temperatures kTi = O.50lo'5i keV and 
fcTa = O.Sto'.i keV, radii Ri = 0.61° ? kin and fla = 0.12l[5:J!; 
km (for a distance of 3.3 kpc: ICaswell et alj|l975l) . and absorption 
A'h ~ (7I2) >< 10^^ cm~^ (Icr errors). The observed averaged 
flux is ~1.7 X 10~^^ erg cm~^ s~ ^, similar to that me asured with 
XMM-Newton in the August 2005 JDe Lucaet ai]|2006h . 

As can be seen from Fig. [T] IE 1613 showed only moderate 
variability between the many Swift pointings, remaining always 
well below the flux level of the 199 9-2000 outburst (~5 x 10"" 
erg cm~'^ s~^: lGarmire et al.l200d) . Also the Chandra HRC-S and 
ACTS -I observations show a flux of~2x 10"^^ erg cm~^ s"^. We 
note that on a statistical ground we are sensitive only to outbursts 
lasting > 1 month. Although shorter phases of enhanced emission 
cannot be ruled out for IE 1613 based on our data, the historical 
behaviour of the s ource indicates that it s outbursts are likely longer 
than a few months jDe Luca et alj2006r) . Thus our data suggest that 
IE 1613 persisted in a quiescent state during the five years of the 
Swift monitoring. 

Owing to their long time span (from MJD 53 804.500 to 
55 632.581), the Swift A&idi are suitable for studying the 24-ks mod- 
ulation of IE 1613. A folded profile showing a significant modu- 
lation cannot be obtained from the XRT data using the most accu- 
rate period available so far, i.e. the one estimated fr om the August 
imS XMM-Newton observation (P = 24.0±0.1 ks: lDe Luca et all 
120061) . This is not surprising since ~200 days separate the XMM- 
Newton observation from the start of the Swift monitoring, and the 
XMM-Newton period uncertainty implies a phase uncertainty of 
half a cycle after only ~30 days. 

So, as a starting point, we computed a fast-Fourier-transform 
power spectrum using all the Swift data at the highest resolution 
allowed by the PC mode (bin time 2.5073 s). Given the approx- 
imate knowledge of the source period, a 'blind' search is not, in 
principle, necessary, but we did this to have a clear picture of the 
Swift time-series: considerable noise can be expected to result from 
both source flux variations and the Swift uneven sampling of the 
light curve, and a search restricted around the XMM-Newton pe- 
riod would have involved the risk of selecting a spurious signal, 
in the case the true signal was embedded in a high level of non- 
white noise. As expected, significant noise is present, but a very 
prominent peak at 24031(2) s (the quoted uncertainty indicates the 



Fourier period resolution) stan ds out well above the noise level, 
with a Leahy-normalised power dLeahv et alJl983r) of 425 (Fig.[3](. 
While the non-white noise does not affect the frequency of a real 
signal, it alter s the statistical proper ties of a time-series; following 
the recipes of llsrael & Stellal ( 11996!) we estimate that, after taking 
into account the number of frequencies searched, the probability 
of having a signal this strong by chance coincidence is lower than 
3 X 10" (that is a detection at a higher than 5.9cr confidence level). 
Moreover, the signal is consistent with the periodicity measured by 
XMM-Newton. We also note that no other periodicity with signifi- 
cance higher than 3cr was found up to ~1 year. 

We folded the Swift data, as well as those of the observations 
in Table [T] on the period P — 24 031 s. We obtained very sig- 
nificant pulse profiles; those of the data taken during the quies- 
cent state of IE 1613 {Swift, XMM-Newton, Chandra/RRC-S and 
ACIS-I) are single-peaked (and well-modelled by two or three sine 
functions with the periods fixed at the fundamental period and 
higher harmonics, with phases and amplitudes free to vary), while 
that from the Cluindra/ ACIS-S observation shows two asymmetric 
peaks per cycle, both exhibiting two sub-peaks. In order to obtain 
a refined ephemeris, we studied the phase evolution through the 
epoch-fo lded data by means of an iterative phase-fitting technique 
(see e.g. iDall'Osso et al.ll2003h . Given the variability of the pulse 
shape, we did not make use of a pulse template to cross-correlate 
with, but we inferred the phase of the modulation by fitting each in- 
dividual folded profile with the fundamental plus three higher har- 
monics[j 

As a first step, we fit the phases of fundamental harmonic 
obtained from the Swift data divided into four segments of ap- 
proximately equal length and from the XMM-Newton and Chandra 
HRC-S and ACIS-I observations. At this stage, the Chandra/ ACIS- 
S observation was left out, because of its very different pulse pro- 
file. The time-evolution of the phase can be followed unambigu- 
ously throughout all the data and described with a linear relation of 
the form (/) — (/!)o + 27r(f — io)/P. We assumed the start of the 5w//f 
monitoring, MJD 55 804.0, as the reference epoch to and the fit 
(xl = 0.88 for 5 dof) gives Pa = 24 030.42(2) s (Icr uncertainty, 
valid over the range MJD 53 605-55 632). We designate this rota- 
tional ephemeris 'solution A'. A quadratic term — 7r(i — to)^P/P^, 
which would reflect the presence of a period derivative (P), is 



•^ This is the minimum number of harmonics necessary to properly fit 
all the folded profiles, including the multi-peaked Chandra! ACiSS one. 
The timing analysis described in the following was performed also with a 
slightly different approach: the phases of the individual folded profiles were 
derived by a fit with a vaiiable number of harmonics, determined on a case- 
by-case basis by requesting that the addition of a further (higher) harmonic 
is not statistically significant (by means of a Fisher-test) with respect to the 
null hypothesis; the results are essentially identical. 
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Phose 

Figure 4. 32-bin epoch-folded pulse profiles of IE 1613 obtained from dif- 
ferent instruments (as indicated in each panel). The blue line distinguishes 
the Chandra/AClS-S data, which were used in deriving solution B but not 
for solution A (see Section[3]for details). In red we over-plotted the XMM- 
Newton/EPlC data before the filtering for proton flares. The fundamental 
harmonic of the pulse profile is shown in green in each panel (in the EPIC 
panel it refers to the filtered data). 



not required. This implies an upper limit on the period derivative 
of IE 1613 of |Pa| < 3.3 x 10"^ s s~^ (3a confidence level). 
In Fig.|4]we show the epoch-folded pulse profiles (including the 
Chandra/AClS-S one) obtained using this solution. We note that 
the XMM-Newton observation was affected by rather intense proton 
flares. The removal of the intervals of flaring background, because 
of the few cycles contained in the observation, significantly affects 
the pulse profiles. As a check of the robustness of our results, we 
repeated the timing analysis using the unfiltered EPIC data and we 
obtained virtually identical results (both filtered and unfiltered pro- 
files are plotted in Fig.|4ll. 

Using solution A, we are able to predict for IE 1613 the phase 
of the fundamental harmonic at the epoch of the Chandra! ACiS- 
S observation within ±0.03 cycles (at 3(t). The phase of the fun- 
damental harmonic measured in the Chandra/ ACIS-S data nicely 
dovetails with the predicted value. Thus we derived a new coherent 
timing solution, which we denote with 'B', including the Chan- 
dra/ ACIS-S data and therefore valid over the range MJD 52 336- 
55 632. Again, the phases of the fundamental harmonic can be fit 
with a linear relation (xt = 0.77 for 6 dof) which yields Pb = 



24030.42(2) s (la uncertainty; epoch MJD 55 804.0). While the 
best-fitting period is equal to that of solution A, the limit on the 
period derivative is slightly more constraining: \Pb\ < 1.6 x 10~® 
s s"'^ (3(7 confidence level). 



4 DISCUSSION 

We presented the analysis of the first five years (2006 April-2011 
April) of the Swift/XKF monitoring of the enigmatic X-ray source 
IE 1613 at the centre of the SNR RCW 103. During this time span, 
the source remained in a quiescent state, with an average observed 
1-10 keV flux of ~1.7 x 10"^^ erg cm"^ s"\ 20-30 times lower 
than the historical maximu m attained in the 19 99-2000 outburst 
(~5 x 10"^^ erg cm"^ s'^i lOarmire et alj200(ih . The timing study 
of the Swift data yielded an accurate measure of the modulation 
period of IE 1613 which allowed us to phase-connect the XRT 
data with archival Chandra and XMM-Newton observations. We 
derived two timing solutions, labelled A and B, both consistent 
with the same constant period Pa = Pb = 24 030.42(2) s but 
with somewhat different upper limits on the period derivative. So- 
lution A (JPaI < 3.3 X 10"^ ss"^ MJD 53 605-55 632) is based 
on the single-peaked pulse profiles observed while IE 1613 was in 



quiescence. Solution B (|Pb| < 1.6 x 10"® s s"^ MJD 52336- 
55 632) is a natural extension of solution A including the multi- 
peaked Chandra/ ACIS-S profile obt ained on 2002 Marc h 03, when 
IE 1613 was in a rather bright state JSanwal et alj|2002n . 

The phase-coherent timing technique employed in our timing 
study of IE 1613 closely parallels the well-tested (virtually all sys- 
tematics are under control) procedures used for X-ray bright spin- 
powered pulsars. For such objects it is implicitly assumed that the 
pulse profile does not intrinsically change in time, so that, when a 
Fourier decomposition of the pulse is introduced, neither the fun- 
damental nor the higher harmonics evolve. In such a picture all the 
variability is due to Poisson noise and the phase evolution can be 
tracked by matching the pulse profile at any given epoch always 
with the same template. 

In the case of IE 1613 there are clear indications that the pulse 
shape changes in time so that the previous assumptions are not 
valid. This forced us to abandon the standard template-matching 
analysis and follow instead the phase of the fundamental (first) har- 
monic|jThe main justification for such an approach is that usually 
(moderate) pulse shape changes are associated with higher harmon- 
ics while the fundamental is stable. Even if this seems to be the case 
for IE 1613, given the stability of the phase of the fundamental har- 
monic over sever al years, we stress th at it does not have to hold in 
general (see, e.g.. lHartman et al.ll2008l) . 

In the following we discuss the implications of our newly de- 
rived upper limit \P\ < 1.6 x 10~® s s^^ on the models that have 
been proposed so far for IE 1613. We remark that all the ensuing 
considerations are based on the assumption (discussed above) that 
the fundamental harmonic is a good tracer of the timing behaviour 
of the source. 

Although the nature of IE 1613 is still an open is- 
sue, most interpretations favour the neutron star scena rio, 
in w hich the star is either isolated JDe Luca et al.l 120061 : iLiI 
l2007h or in a bina r y system with a low-mass companion 
jPizzolatoetai] l2008l: iBhadkamkar & GhoshI l2009h . While the 



^ The large phase uncertainties do not allow a similar study of the higher 
harmonics (which are not even always detectable). 
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model bv lBhadkamkar & GhoshI ( 120091) is based on a fast-spinning, 
moderately magnetised neutron star in a 6.67 hr eccentric orbit, in 
all the other cases an ultra-magnetised neutrons star (B ~ 10^^ G) 
is required to explain the very long period of IE 1613 and the ob- 
served X-ray peri odicity is related to th e star spin period (in the 
binary scenario of IPizzolato et alJuOOa the latter may or may not 
coincide with the o rbital period). Actually, as already noted by 
lOe Luca et al.l ( |2006h . magneto-dipolar braking alone is not enough 
to explain the present value of P even invoking a magnetar, and 
spin-down by the interaction of the star magnetosphere with a 
(residual) disc is also necessary for an isolated object. Spin-down to 
P ~ 6.67 hr in a time ~2 ky r can be achieved if th e initial period 
is peculiarly long (>300 ms, lOe Luca et al. 20061: magnetars are 
in fac t believed to be born with ms periods, [Thompson & Duncarj 
1 19930 . U sing different assum ption s on th e termination radius of the 
disc (see lEksi & AlpailbOOSh .lLiI ilOOH ) showed that the same re- 
sult can be recovered also for more conventional values of the initial 
period (~10 ms). 

The current spin-down rate expected in the binary magnetar 
model is very small since the equilibrium period is reached well in 
advance of ~ 1000 yr, unless the synchronisation time is very s hort 
(~10 yr) and there is no mass transfer in the system (Pizzolato e t al.l 
uOOa) . If IE 1613 is an object of this kind our upper limit on P is 
completely non-constraining for the model. 

On the other hand, if IE 1613 is an isolated magnetar sur- 
rounded by a fossil disc it must be a quite rare system. According 
to the Monte Carlo simulations of Ilj 1 120071) , the fraction of objects 
of this type with periods >100 s at an age of 2.5 kyr is only ~1%. 
The large majority of stars are much faster rotators (P ~ 1-10 
s) in the ejector stage (typically identified with SGRs/AXPs) while 
ultra-slow systems are in the propeller/accretor phase[j 

The spin-down rate of a neutron star which inter- 
acts with a fossil disc in the propeller stage is given by 
P - -MR^nK{R^n) - 2-k/P]P'^/{-kI) where i?,„ = 
0.5[B''7i^7(8G'MM2)]i/^ is the Alfven radius, M - 1.4Mq 
and R ^ 10^ cm are the star mass and radius (/ ~ 10*^ g cm'^ 
is the moment of inertia), B is the surface magnetic field, SIk is 
the Keplerian angular veloc ity and M oc t^^'^^ is the mass loss 
rate from the disc JLil 120071 . and references therein)|j If IE 1613 
is currently in the propeller phase, P can be derived from the pre- 
vious expressions with P ~ 24 ks and i ~ 2.5 kyr, as a func- 
tion of the disc initial mass Md{0) and of B. Results are shown 
in Fig. [5] (lower panel) for 10"^ < Md{0)/Mq < 10"^ and 
3.2 X 10^^ G < B < 10^*^ G; the Alfven radius as a function 
of Md (0) and B is plotted in the upper panel. Comparison of Rin 
with the light cylinder radius, Rlc = cP/{2n), and the co-rotation 
radius. Re = [GMP'^ /{Atv'^)]^^'^ shows that the star is currently 
in the propeller phase, characterised by Re < Rin < Rlc, only if 
B > 10^'' G. Lower values of the field or large enough (depending 
on B) initial disc mass would result in IE 1613 being in the accre- 
tor stage. The resulting values of P are orders of magnitude above 
the upper limit we r eported (as a reference, the case discussed in 
lOe Luca et al.ll2006l has P ~ 10^^ s s~^), with the only exception 
of an extremely narrow range of Md (0) for each B. Although not 
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Figure 5. Upper panel: the Alfven radiu.s as a function of the initial disc 
mass AloiO) for different values of the magnetic field in the range 13.5 55 
log B (G) s; 16 (the field increases as the curve moves upwards); the light 
cylinder radius R^c ^id th^ corotation radius Re are also shown. Lower 
panel: the period derivative as a function of Mj^ (0) for those values of 
B for which the source is cun'ently in the propeller stage; the dashed line 
marks our derived upper limit on P. Here the period and age have been 
fixed to those inferred for IE 1613 (see text). 



all the values of AId (0) and B we considered will produce the cor- 
rect period at an age of 2.5 kyr (this depends also on both the initial 
period a nd the angle between the magnetic and spin axes, see again 
lLill2007n . our conclusion is that if IE 1613 is an isolated magnetar, 
and hence it must be either a propeller or an accretor at present, 
the typical large expected values of P 2> 3 x 10~® s s ~^ makes 
the propeller option rather unlikely. Assessing the reliability of the 
accretor scenario requires a more thorough analysis. Here we just 
note that if the same expression of the torque remains valid during 
the accretion phase, the spin-up rate is also largely in excess of our 
upper limit on P. 



^ Although previous figures are model-dependent the conclusion that pro- 
pellers/accretors are a tiny minority appears to be robust. 
^ We remark that no complete theory of the intera ction of a disc with a 
magnetised neutron star exists. Following iLiI (2007) we adopted the 'effi- 
cient' form of the propeller torque (see e.g. [Francischelli & Wiieral2002l 
for a comparison of different expressions of the torques). 
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